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ABSTRACT 

We study a nonlinear development of radial pulsation instability to a resonant excitation 
of nonradial modes. Our theory covers the cases of axisymmetric (m = 0) modes as well as 
(m, — m) pairs. Adopting a simplified treatment of the radial and nonradial mode coupling 
we find that the asymptotic state is a pulsation with constant amplitudes and we evaluate the 
relative amplitude of the nonradial component. 

Observable consequence of the m = mode excitation is a small period change and a more 
significant amplitude change, especially in the case of a dipole mode (1 = 1). Such a mode has 
a fairly large excitation probability. 

Significant amplitude and phase modulation is predicted in the case of excitation of a 
m = ±l pair. We suggest that this may explain Blazhko-type modulation in RR Lyr stars. 
If this model is correct, the modulation period is determined by the rotation rate and the 
Brunt- Vaisala frequency in the deepest part of the radiative interior. 

Stars: oscillations - Stars: variables: RR Lyr 



1 Introduction 

Among a variety of types of pulsating stars, RR Lyr variables seem to be the 
most important ones. They are significant in various fields of astronomy: dis- 
tance measurements, evolution of Population II and low mass stars, structure 
of Galaxy and, obviously, stellar pulsation theory. Thus, it is important to 
understand the nature of all phenomena connected with their pulsations. One 
can distinguish a few subclasses of RR Lyr stars: RRab stars pulsating in the 
fundamental radial mode with large amplitude; RRc stars - the first overtone 
pulsators with small amplitude and sinusoidal light curves; RRd stars that pul- 
sate in both fundamental and first overtone modes. 

The most intriguing phenomenon that has been waiting for theoretical expla- 
nation for almost a century is the Blazhko effect. It is characterized by periodic 
modulation of oscillation amplitude and phase on a longer time scale. The effect 
was discovered by Blazhko (1907). Basic properties of observed Blazhko-type 
modulations are summarized by Szeidl (1988) and Kovacs (1993, 2000). Accord- 
ing to these authors, about 20-30% of all RRab stars exhibit this phenomenon. 
The modulation is often quite strong, pulsation amplitude may change by a fac- 
tor of up to 3 during the Blazhko cycle. Fourier analysis of light curves with the 
Blazhko effect shows two peaks, symmetrically spaced on both sides of the main 
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peak corresponding to the basic oscillation frequency, as well as similar peaks 
around harmonics of the main frequency. This feature of Fourier spectrum is 
by some treated as a definition of the phenomenon. However, we should stress 
that the Fourier analyses are not available for most of RRab stars reported to 
show the Blazhko effect. 

Whether RRc stars exhibit this type of modulations was a matter of spec- 
ulations for years. Recently such variables were discovered in the LMC by the 
MACHO project (see Kurtz et al. 2000, Alcock et al. 2000). However, incidence 
of the Blazhko phenomenon among RRc stars is much lower than in the RRab 
case and is estimated at about 2%. 

Early explanations of the Blazhko effect postulated the existence of strong 
magnetic field with the axis inclined to the axis of rotation. This field was 
thought to cause a nonspherical distortion of the pulsation. Stellar rotation 
led to the amplitude modulation seen by observers. A mathematical model of 
magnetically modified pulsation was developed by Shibahashi (1995) (see also 
Shibahashi 2000). There are many basic problems connected with this model 
(see Kovacs 1993, 2000) and we will not deal with it in this paper. 

Alternative interpretations invoke nonlinear interactions between oscilla- 
tion modes. Two approaches have been adopted: numerical simulations of 
fully-nonlinear pulsation and an approximate formalism of amplitude equations 
(AEs). The latter is the only available formalism for nonradial oscillations. It 
proved particularly useful to study effects of resonances (see e.g., Dzicmbowski 
1982, Buchler and Goupil 1984, Van Hoolst 1992, 1994). The AE-formalism 
finds its applications in various fields of physics. Our application bears most 
resemblance with work of Wersinger et al. (1980) who studied nonlinear interac- 
tion of plasma waves in the case of parametric resonance. These authors showed 
that solutions of AEs may be stationary (fixed-point), periodically modulated 
(periodic limit cycle) as a result of Hopf-bifurcation, or even chaotic. In the 
case of pulsating stars parametric resonance is likely responsible for limiting of 
oscillation amplitudes in <5-Sct stars (Dziembowski and Krolikowska 1985). 

Moskalik (1986) proposed 2:1 resonance between the fundamental and the 
third overtone modes as a possible explanation of the Blazhko effect which would 
be a manifestation of the periodic limit cycle solution of the corresponding AEs. 
However, the values of parameters necessary to obtain modulated solutions seem 
to be outside the range of RR Lyr models (Kovacs 1993 and references therein). 
Moreover, there is no sign of modulation in hydrodynamic simulations (Kovacs 
and Buchler 1988). 

Nowadays, the most promising hypotheses involve excitation of nonradial 
modes. Properties of such modes in giant pulsators (RR Lyr, Cepheids) were 
studied already over twenty years ago (Dziembowski 1977a, Osaki 1977) but 
until the early 90s these papers were largely out of interest of astronomers. 

Van Hoolst (1992) derived AEs for resonance between two modes with nearly 
equal frequencies (the 1:1 resonance) in the adiabatic approximation. He cal- 
culated polytropic stellar model and found a nonradial mode with frequency 
close to that of the second radial overtone. For these two modes he calculated 
some parameters in AEs. For other parameters he considered wide ranges of 
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values. Depending on the values of parameters he found different types of solu- 
tions which include single mode pulsation, double mode pulsation with constant 
amplitudes, pulsation with periodically or chaotically modulated amplitudes. 

Kovacs (1993) and Van Hoolst and Waelkens (1995) proposed the 1:1 reso- 
nance between the radial fundamental mode and a nonradial mode as a possible 
explanation of the Blazhko effect. This type of resonance is not possible be- 
tween low order radial modes. Van Hoolst et al. (1998, hereafter VDK) showed 
that in realistic RR Lyr models, nonradial mode spectra in the range of the 
fundamental and the first overtone radial modes are very dense. Hence the 
frequency distance between the radial and the nearest nonradial mode is very 
small. VDK studied stability of monomode radial pulsation in a representative 
RR Lyr star model. They found high probability of resonant excitation of the 
nonradial mode. The highest probability occures for modes with / = 1,5,6 and 
1 — 1,4 at fundamental and first overtone modes, respectively. Dzicmbowski and 
Cassisi (1999, hereafter DC) repeated similar calculations for a set of realistic 
RR Lyr models covering the whole range of effective temperatures corresponding 
to instability strip in RR Lyr stars luminosity region and found that conclusions 
of VDK are generally valid. They found, in particular, that the probability of 
resonant excitation of nonradial modes is about 0.5, which is not far from the 
observed incidence rate of the Blazhko effect among RRab stars. According to 
their calculations, the effect should occur among first overtone pulsators nearly 
as often as among RRab stars, which is in a strong contradiction with obser- 
vations. However, they stressed that excitation of a nonradial mode does not 
necessarily imply amplitude modulations. 

Nonlinear consequences of nonradial mode excitation due to the 1:1 reso- 
nance have never been studied in application to RR Lyr stars. The Blazhko 
effect may be a manifestation of periodic limit cycle solution of AEs, but it 
may also be caused by nonradial mode excitation with constant amplitude. The 
latter possibility requires m^O. In this case, the Blazhko period is related to 
but not identical with the rotation period. 

In this paper we extend the formalism of AEs used by VDK to study solutions 
with finite amplitudes of nonradial modes. We keep simple scaling of resonant 
coefficients adopted by DC. 

The plan of the paper is as follows. In Section 2 we outline the basic linear 
properties of nonradial modes in evolutionary models of RR Lyr stars. The 
properties are described in details by VDK and DC. In Section 3 we transform 
AEs of VDK into the form used in further part of this paper, by adopting 
some simplifications, e.g., above mentioned scaling of coefficients. We study 
stationary solutions and their stability in Section 4, as well as time dependent 
solutions and attractors in Section 5. Section 6 is devoted to estimating the 
mode amplitudes and period changes caused by the resonant interaction in all 
DCs models, as well as the observational effect of m = mode presence. In 
Section 7 we explain how the (m,—m) pair excitation may cause the observed 
amplitude and phase modulation. Section 8 includes conclusions and discussion 
of the results. 
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2 Nonradial Modes in RR Lyr Models 

RR Lyr stars are Horizontal Branch helium burning giants, characterized by 
strong mass concentration in the center. VDK showed that in such stars the 
oscillation equations can be solved analytically in the deep interior, using an 
asymptotic approximation. The solution may be matched to the numerical 
solution in the envelope. Applying this procedure one can determine basic 
properties of nonradial modes. The most important is that there is always 
dense spectrum of nonradial modes near the two lowest degree radial mode 
frequencies. For a given I, the mode frequency decreases with increasing n. This 
is so because, owing to very high values of the Brunt- Vaisala frequency in the 
interior, all nonradial modes are in fact (/-modes of high degree. When rotation is 
ignored all modes with the same l,n have the same frequency. Non-zero rotation 
means still denser spectrum for each I. In the envelope nonradial modes closest 
in frequency to a given radial mode have their radial eigenfunctions almost the 
same as the radial mode eigenfunctions. This remains true as long as the mode 
frequency is much higher than the Lamb frequency (L; = c(l + 0.5)/r where c is 
the local sound velocity). 

Radial modes are of course purely acoustic and they propagate only in the 
envelope, where the Brunt- Vaisala frequency is lower than their eigenfrequen- 
cies. The outer layers determine properties of radial modes and the interior 
plays no role. 

From linear nonadiabatic calculations we find time dependence of pertur- 
bations in the form exp(io>j +Kj) where u)j is the eigenfrequency and Kj is the 
linear growth rate of mode numbered by j. The growth rate may be expressed 
as (see e.g.,Unno et al. 1989) 

-- ^- <u 



J MM 

where Ij is the inertia of the mode and Wj is the work integral, i.e., the integral 
over star's volume of nonadiabatic operator acting on the eigenfunction. In the 
case of the radial mode it is enough to integrate over the envelope, because the 
contribution from the deep interior is negligible. In the case of the nonradial 
modes one may calculate separately the integrals over the deep interior (gravity 
waves propagation zone) and over the envelope (acoustic waves propagation 
zone). Moreover, the eigenfrequencies of both modes are nearly the same, and 
one can put uj instead of uji in the formula for the nonradial mode growth rate. 
With above assumptions one obtains 

° 2/ wo' ^ ^ 

W + W g 

Kl = 2I{U = Kl ' + K a- ( 3 ) 

The quantity n g = W g /(2IiWo), which is always negative, measures the damping 
rate in the deep interior. 

DC calculated three sets of evolutionary models of RR Lyr stars (for differ- 
ent masses and metallicities) . They used the procedure described by VDK for 
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these models and calculated properties of the two lowest radial modes and the 
corresponding resonant nonradial modes for a few lowest Vs. In the vicinity of 
the fundamental mode frequency, nonradial modes with I > 7 are damped very 
strongly and they are never excited due to 1:1 resonance with the radial mode. 
The same is true for the modes with I > 6 in the vicinity of the first overtone 
frequency. 

In all numerical calculations in this paper we use data for models from DC. 



3 Amplitude Equations 

The amplitude equations used in this paper are based on those derived by 
Van Hoolst (1994) and used by VDK in their study of stability of radial pul- 
sation in RR Lyr stars. Here we first derive our simplified form for the case of 
resonant coupling between a radial (/ = 0) and an axisymmetric nonradial mode 
(I > 0, m = 0). Then we will derive corresponding equations for the case of cou- 
pling to a nonradial mode pair of the same degrees I and opposite azimuthal 
degrees, m and — m. We will see that there are no essential differences between 
these two cases. Some details of the derivation are deferred to Appendix A. 

We introduce complex pulsation amplitudes a associated with an individ- 
ual oscillation mode (j), which are defined by the following expression for the 
displacement 

Sr 1 

— =&(r,t) = -aje^wj^ + cc, (4) 
' o ^ 

where ro is the Lagrangian radial coordinate, Wj is vector eigenfunction, and 
LOj is eigenfrequency of the mode. We adopt normalization 

Wj, r =Y t ^ mj \/4Tr (5) 

at the surface, which implies that \a\ is the rms amplitude of the surface radius. 



3.1 The m = Case 

The time dependence of amplitudes aj = aj exp(iu) jt) is described by AEs. For 
the 1:1 resonance between the radial mode and axisymmetric nonradial mode 
the relevant AEs are (VDK) 

^ = (K + iwo)ao + 5oao|ao| 2 + 5oa |a;| 2 + 

+R 1 Q a*a 2 + R 2 a l \a l \ 2 , (6) 

^ = (K l + iLu l )a l + Sl > a l \a \ 2 + S l l a l \ai\ 2 + 

+Rl^a 2 + Rf\ ai \ 2 ao + Rfa 2 a*, (7) 

where Sj, Rj denote saturation and resonant coupling coefficients, respectively. 
For odd I coefficients R 2 , Rf and Rf vanish. In our work we devote special 



6 



attention to the case of I = 1 and thus we neglect all the terms involving these 
coefficients. We believe that application of so simplified system to even Z's will 
not result in large error because these terms are of higher order in the amplitude 
of the nonradial modes. We will see that the amplitudes of I > 1 modes are 
always much smaller than that of the radial mode. We can now introduce real 
amplitudes A, and phases <f>j in the form aj = Aj exp(iujjt + <j>j) and separate 
Eqs. (6) and (7) into the real and the imaginary parts. This yields 

dA 



dt 

dA 
dt 



A a(Po 



dfr 



K +n(s^A 2 +n(s l Q )Af + 

-(3?(i?o) cosr - %( R o) sinr ) A f 



«, + R(S?)4> + H(S/)A? + 
+ (n (Rj ) cos r + S(R} ) sin r) Al 1 Ai , 

Z(S )At + Z(S l )A* + 
+ (5R(J?J)sinr + 9( J Rj)cosr)^ 

Q(S?)At + Q(S l l )Al + 



-(-3t(R})sinT + Q(R})cosr)A 



where we introduced the relative phase, 

r = 2(Awt + <j>i-<f> ), 



(8) 



(9) 



(10) 



(11) 



(12) 



and the detuning parameter, Au> = wi — loq. The time derivatives of the phases 
4>j give nonlinear corrections to the linear mode frequencies. Terms with satura- 
tion coefficients in Eqs. (10) and (11) give nonresonant corrections and may be 
neglected because they move the whole spectrum of nonradial modes without 
significant change of frequency distance between consecutive nonradial modes. 
Moreover, we assume that the resonant coupling coefficients may be regarded 
purely imaginary as in the adiabatic case (Van Hoolst 1992). Although pulsa- 
tion is never strictly adiabatic, we follow here a suggestion of VDK who argued 
that real parts of these coefficients are much smaller than imaginary parts. 

All coupling coefficients are proportional to the complicated integrals con- 
taining products of four eigenfunctions for the modes involved. These integrals 
can be separated into integrals over radial and angular coordinates. The latter 
are easy to calculate analytically, because they contain well known spherical 
harmonics. The former are not possible to be determined analytically. Never- 
theless, some useful relations between coefficients can be derived. This is due 
to the fact that the radial eigenfunctions of both radial and nonradial modes 
are almost the same in the outer layers and integration of products of different 
eigenfunctions gives the same result. In Appendix A.l it is found, in particular, 
that 



= 25°, 



(13) 
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where Hi= j P l (x)dx and P ; is normalized Legendre polynomial. The values 

-l 

of Hi for a few low I's are given in Table 1. 

Table 1 
Numerical values of Hi 



I 


Hi 


1 


0.9 


2 


1.071 


3 


1.180 


4 


1.260 


5 


1.322 


6 


1.373 



Like the growth rates (see Eqs. 2 and 3), all the coupling coefficients are 
proportional to the inverse of mode inertia. Thus, they satisfy the relations 



R(sS)_»(sf) = _ rt 

— CI, 



K 



(14) 



(15) 



Let us note that a determines the amplitude of single-mode radial pulsation. 

When we introduce the dimcnsionless time T=n t, the dimensionless fre- 
quency distance Aa = Aoj/kq, and the following quantities 



7 = 



Kin 



Kl,0 
K ' 



(16) 



we can, with the help of Eqs. (2), (3), (14), and (15), rewrite Eqs. (8)— (11) in 
the form 

dA 

dr 
dAj_ 

dr 
At 

We neglected the real parts of the resonant coefficients and the imaginary parts 
of the saturation coefficients, as was explained above. From the last two equa- 
tions and Eq. (12), with an additional assumption of non-zero values of both 
amplitudes, we obtain 

dr ^ ^)=2{A„ + RcosT(lAl-A>)). (21) 





-RsiaTAf] A , 


(17) 


l[l- 1 -2a{Al- 


^HiA^+RsinTAl] A u 


(18) 


RcosTAfA , 




(19) 


IRcosTAlAi. 




(20) 



dr V dr 
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Eqs. (17), (18), and (21) form a complete set of equations and we will use 
them to study double- mode solutions. The equations are invariant under the 
transformation 

(Aa,r)-(-Aa,7r-r). (22) 

Thus, we will restrict ourselves to Act > 0. 
3.2 The m^O Case 

Resonant interaction between radial and nonradial modes with nearly equal 
frequencies is only possible when the nonradial mode has azimuthal number m 
equal zero. In other cases, all resonant coupling coefficients vanish, as was shown 
by Van Hoolst (1994). Nevertheless, modes with m^O can interact with a ra- 
dial mode due to the resonance between the radial mode and a pair of nonradial 
modes of the same value of I and opposite values of m. Analogical situation of 
resonance between modes of rotationally split triplet I = 1 is described by Buch- 
ler et al. (1995). The difference is that in our case the m = mode is replaced by 
the radial mode. These authors call this kind of resonance the 1:1:1 resonance 
but we think a more appropriate name is the 2:1+1 resonance. The reason is 
that the three frequencies satisfy the resonant condition 2w ~w m +w- m - 

The AEs for the above-mentioned type of resonance have the form given 
e.g., by Buchler et al. (1995). The relations between coupling coefficients similar 
to those in the m = case are derived in Appendix A. 2. The relative phase T 
in this case is given by 

r = 2Aujt + (/>_ + 0+ - 20o (23) 

where Aw = (cj_ + w+)/2 — luq. Such a definition of T and Au> is consistent with 
the m = case. In further studies we will assume that rotation is slow and except 
for the linear eigenfrequencies it does not change the dynamical properties of the 
nonradial modes. The radial eigenfunctions of these modes are m-independent. 
Thus, we have a full symmetry between coefficients in equations for amplitudes 
and phases of both nonradial modes. This means, in particular, that 

K-=K+ = K U S-=S+ = Si, R-=R+ = Ri (24) 

where quantities with subscript I are the same as in the m = case. We stress 
that in the present case there is only one non- vanishing resonant coefficient also 
for even Vs. Similarly as in the m = case we introduce dimensionless quantities 
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and we obtain 

dAo 
dr 



[l - a(Al + 2A 2 _ + 2A\) -2RsinTA-A + ] A , (25) 
/ [1 - 7- 2a (A 2 + HJ n (A 2 _ + 2A 2 + ))]A_ + IRsmTA 2 A+, (26) 
/ [1 - 7- 2a {Al + H\ n (2A 2 _ + A\))\ A+ + IRsmTA 2 A- , (27) 



dr 

r r-i - .).. / |2 rr»'/.) |2 ( 2 ,\" t rn ;,.r (2 

dr 
d^o 
dr 
i(j>- 
dr 
d0+ 
dr 



Ao-p = 2i?cosr^A+A , (28) 
dr 

AJ^- = IRcosTA 2 A+, (29) 
dr 

A+^p- = IRcosTA 2 A_, (30) 



1 

where H™ = J Pj i m (x)dx, P^ m is normalized associated Legendre function and 
Eqs. (23) and (28)-(30) allow us to write equation for the phase T 



^ = 2Aa + RcosT 
dr 



(31) 



Eqs. (25)-(27) together with Eq. (31) form a complete system of equations which 
allow us to determine the amplitudes of the modes participating in oscillations. 

When we multiply Eq. (26) by A_, Eq. (27) by A + , and subtract them, we 
obtain 

-j^ (AL - A\) = 21 [1 - 7- 2aA 2 - 2aH\ n {A 2 _ + A\)} (AL - A\) . (32) 

It will be seen further that a typical value of Aq for stationary solutions is about 
1/a (exact equality is for single- mode solution). This means negative value of 
the square bracket in Eq. (32) and exponential decrease of the difference of the 
two amplitudes. Thus, we limit ourselves to the A_ =A + solution. 

Let us now focus on the case of modes with 1 = 1, m = ±l. Then we have 
Hi =3/5. When we introduce a pseudo amplitude A\ 

A-=A + ^^=. (33) 

then Eqs. (25)-(27) and (31) are identical with Eqs. (17), (18), and (21), where 
.Hi = 9/10. Thus, we may consider the pair of nonradial modes 1 = 1, m = ±l 
as an equivalent of the mode 1 = 1, m = 0. This fact has a simple physical and 
geometrical explanation: spherical harmonic Y^o is equivalent to the sum of 
spherical harmonics Yi,-i + 3'i,i divided by \pl in the system of coordinates 
with the z-axis inclined by ir/2. The pair of modes 1 = 1, m = ±l excited with 
the same amplitudes form a distortion of the star that has the same shape 
as the distortion associated with the mode 1 = 1, m = but with the axis of 
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symmetry lying in the plane of equator. Since all physical mode characteristics 
are m-independent such a pair interacts with the radial mode in the same way 
as the 1 = 1, m = mode. Therefore, the amplitude Ai determined as a solution 
of Eqs. (17), (18), and (21) gives also the amplitudes of the nonradial modes 
A ± =A t /V2. 

We treat the interaction of the radial mode with the 1 = 1, m = mode and 
with the 1 = 1, m = ±l pair independently. Strictly speaking this is justified 
only for the instability analysis. The small quadratic effects in the rotation rate 
cause that (wi,i + u;i j _i)/2 7^u;i i o and these two cases are indeed independent. 
However, at finite amplitudes of the nonradial modes the interaction between 
the modes within the triplet is likely to be important. It is possible that the 
whole triplet will be exited but this should not have large consequences for 
amplitude prediction. 

Similarity of the 1:1 and the 2:1+1 resonances, involving the (m, — m) dou- 
blet, extends to I > 1. The only difference concerns the Hi coefficient (see Eqs. 26 
and 27). 

It is also possible for a pair of nonradial modes belonging to different multi- 
plets (i.e., of different radial orders) to be excited. Then the the amplitudes of 
the two modes cannot be assumed equal. The difference is mainly due to the 
difference in mode inertia, which is strongly frequency sensitive. We define I± 
similarly as in Eq. (16), that is we set 

/ ± = T^, 04) 
and, ignoring the difference in the coefficients at H™, we get 

(35) 

If this ratio is not too different from unity then, upon introducing 




7= v ^O+, Af = 2A + A_, (36) 

we reduce equations describing such interaction to Eqs. (17), (18), and (21) 
which describe the two-mode resonance. 



4 Stationary Solutions 

Stationary, or fixed-point solutions of the amplitude equations are obtained 
from Eqs. (17), (18), and (21) where we set dA,/dr = 0. Then the nonlinear 
frequencies are given by u>j -+d<j>j/dT, and, like the amplitudes, are constant. 

4.1 Single-Mode Solutions 

Eqs. (17), (18) admit single mode solutions 

(A° ) 2 = -, (A°) 2 = 0, (37) 
a 
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and 



(A o r = 0, (A l ) 2 = ^. (38) 

It is of some interest to investigate stability of these solutions in the absence 
of resonances. Perturbation of Eqs. (17) and (18) with R=0 leads to a simple 
eigenvalue problem for the growth rate /3. We find for the radial mode solution 
j3 = — 2 and — (I + 7). Thus, this is a stable solution. For the nonradial mode 
solution we find /?=— 1 — 37 and 1 — (1 — 7) /Hi. The latter eigenvalue is less 
than zero, if 7<1 — Hi. Since 7>0, this may happen only for 1 = 1 because 
as we may see in Table 1 Hi > 1 for I > 1. This means that in principle we 
could expect that some of RR Lyr stars are I = 1 pulsators. However, we find 
it unlikely because even if 7 is very small the growth of the I = 1 mode is much 
slower than that of 1 = 0. 

A survey of stability of radial pulsation to excitation of resonant nonradial 
modes in RR Lyr star models was given by DC. However, there was a need to 
repeat the calculations because of an error in the adopted form of the stability 
criterion. The correct criterion derived by VDK is 

Ac 2 >\R}\ 2 (A») 2 -g 2 (39) 

where 

g = Ki + %(S?)(A° ) 2 (40) 

is the nonlinear growth rate of the nonradial mode, and other symbols are the 
same as in Eq. (9). In our dimensionless notation Eq. (39) has the form 



Aa 2 >I 2 



■2 



R 

TT-(l+7) 



(41) 



where we used Eq. (37) and the nonlinear growth rate of the nonradial mode 
is easily obtained from Eq. (18). In statistical approach, the above criterion 
implies that the probability of instability of radial-mode fixed-point solution is 
given by 

p _ (42) 

A(T ma x 

where Acr max is the maximum absolute value of the detuning parameter, which 
in the case of dense spectra of nonradial modes in RR Lyr stars is given by a half 
of the frequency distance between consecutive nonradial modes, (07. „ — <7;. n +i)/2. 

DC calculated probabilities of excitation of nonradial modes in the three 
sets of RR Lyr models. However, they used erroneously — J7 (in our notation) 
as the nonlinear damping rate of the nonradial mode. From our AEs we get 
the damping rate —1(1 + 7), which implies lower probabilities of excitation of 
nonradial modes, in particular in those cases where 7 is very small and a is 
comparable to R. This is what happens in the first overtone vicinity. 

We calculated these probabilities once again for the DC RR Lyr models using 
the i?,7,7, Acr max parameters inferred from the data given in DC. We adopted 
a = 200 for the fundamental mode pulsation and 1200 for the first overtone, 
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Fig. 1. Probabilities of excitation of nonradial modes with 1 = 1,2 in the models of DC. Solid 
and empty symbols denote modes in first overtone and in vicinity of fundamental mode, 
respectively. Encircled symbols indicate I = 2 (only for fundamental mode vicinity) . Symbols 
are the same as in DC and refer to evolutionary tracks of different initial parameters. Squares 
and circles denote models characterized by ZAHB composition Z = 0.001, Yhb = 0.243 and 
masses 0.65 and 0.67 Mq, respectively. These quantities for models denoted by triangles are 
Z = 0.0002, Y HB =0.24 and 0.74 M . 



which, according to Eq. (37), correspond to Aq — 0.0707 and A^ = 0.0289, re- 
spectively. These values may be regarded typical for RRab and RRc stars, 
respectively, and we will use them for all numerical calculations in this paper. 
Our calculations confirm DCs result that the most probably unstable are the 
1 = 1,5,6 modes in the fundamental mode vicinity and the 1 = 1,4 modes in the 
first overtone vicinity. However, in the former vicinity we now find relatively 
higher probability of excitation for the / = 2 modes. The probabilities are pre- 
sented in Figs. 1 and 2. 

4.2 Double-Mode Solutions 

When all amplitudes are constant and non-zero, we also have constant T, as 
can be seen in Eqs. (17), (18), and (21) or in Eqs. (25)-(27) and (31). Constant 
T implies the nonlinear frequencies to fulfill exactly the resonance condition, 
i.e., phase-lock phenomenon. For the 1:1 resonance it means the equality of 
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Fig. 2. The same as in Fig. 1, but for the modes with Z = 4 (overtone), 1 = 5 (fundamental) 
and 1 = 6 (fundamental, encircled). 



the two frequencies, while for the 2:1+1 resonance it means that the three 
frequencies are equidistant with the radial mode frequency being in the center 
of the triplet (see. Eqs. 12 and 23, respectively). Buchlcr et al. (1997) give 
general discussion of the phase-lock for other resonances. 

Using Eqs. (17), (18), and (21) we can express the amplitudes and the de- 
tuning parameter in terms of the relative phase as follows 

M dm,2 _ ^sinr dm (l - 7 ) - 2a(H l - 1 + 7) 
[ ° ' 2*2(2- If,) -(£sinr dm )2 ' 

M dmx2 _ a(l + 7)-flsinr dm 

Act = J?cosr dm [(Af m ) 2 -/(4 m ) 2 ]. (45) 

In general it is impossible to obtain an explicit expression for T dm as a function 
of Act. Thus, fixed-point solutions have to be studied numerically. However, 
one important property of the solutions can be obtained analytically. Let us 
consider transition between single- and double-mode fixed-point, i.e., the vicin- 
ity of Af m = 0. Then i?sinr = a(l + 7 ) and A% m = A% (see Eqs. 43, 44, and 37). 
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With the help of these expressions we get from Eq. (45) 

|Act| = A<^ — ^l-_(l + 7 )2 = 7 ^__(i +7) 2. (46) 

Note that it is the same critical value of Aa as given by Eq. (41). It means 
that the nonradial mode appears with zero amplitude at the onset of instability. 
Further, it may be shown analytically that for |Act| > Aad there are no double- 
mode solutions for which Ai/A <^l. The analysis is not reproduced here but 
the property will be seen in our numerical examples. 

The value of Aad is critical for the consequences of the resonances consid- 
ered here. For \Aa\> Aad the single-mode fixed-point solution is stable. For 
lAcr^AfTd it is unstable, but a double-mode fixed-point occurs. It will be 
shown that this double-mode fixed-point solution is almost always stable. 



4.3 Stability of Double-Mode Fixed-Point Solution 

Let us note that Eqs. (17), (18), (21) can be presented in the form 



dAf 
dr 
cL4] 

dr 
dsinr 

dr 



2(1- aA% - 2aAf - RsinTAf ) A 2 , 
^ = 2l(l-~f-2aA 2 -2H ia A 2 + RsmTA 2 ) A 



= 2cosr[Acr + i?cosr(J^-Af)] . 



(47) 
(48) 
(49) 



Let us now consider a small perturbation of the double-mode fixed-point solu- 
tion, 

Al = (A d ™f + 5Al A 2 = {A^f + 5Al (50) 



sin T = sin Tdm + <5 sin T , cosr = cosrd m + £cosr. 



(51) 



Now we insert the above quantities to Eqs. (47)-(49) and neglect terms of second 
and higher orders with respect to small perturbations. In this way we obtain a 
system of linear differential equations which can be written in a matrix form 



dX 

dr 



= AX 



where 



and 



A = 2 



X = 



5AI 
SAf 
dsinT 



(52) 



(53) 



-a A 2 



I(-2a + RsmT)Af 
7i?(l-sin 2 r) 



-(2a + RsmT)Al 

-2IHiaAf 
-i?(l-sin 2 r) 



-RA 2 A 2 
IRA 2 A 2 



-RsmT(IA 2 -A 2 ) 



(54) 
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We made use of cos 2 T = 1 — sin 2 T , cos T S cos r = — sin T 5 sin T relations and omit- 
ted sub- and superscripts dm denoting double-mode fixed-point solution. 

Eq. (52) implies time dependence of the perturbations in the exponential 
form, exp(Ar). There arc three eigenvalues determined by the following cubic 
equation 

A 3 + a 1 \ 2 + a 2 \ + a 3 = (55) 

where 

ai = a(A 2 l + 2IH l A 2 ) + Rsmr{IA 2 a -A 2 ), (56) 
a 2 = Ia 2 AlA 2 l (2H l -A) + aRsmT{IAl-A 2 l ){Al + 2IH l A 2 ) + 

+/i? 2 A 2 A 2 (2-sin 2 r), (57) 
a 3 = IA 2 Af[RsmT(IA 2 -Af)(R 2 -a 2 ('i-2H l )) + 

+R 2 (1 - sin 2 1> (A 2 (l + 21) + 2A 2 (1 + IH t )) ] . (58) 

The stationary solution is stable when all three eigenvalues have negative real 
parts. The condition for the stability is then given by Hurwitz criteria 

ai > 0, (59) 
aia 2 -a 3 > 0, (60) 
a 3 > 0. (61) 

We will check these inequalities numerically to determine stability of each of 
our stationary solutions. 



4.4 Numerical Examples 

In this Section we will show examples of double-mode fixed-point solutions for 
the models from DC and will discuss their stability. The solution is characterized 
by the parameters 7, R, a, I, Act. For the detuning parameter Act we will 
consider the range of values from to A(T max . 

We choose representative examples of nonradial modes which have the high- 
est excitation probability in the vicinity of the fundamental mode and the first 
overtone. 

In Figs. 3-5 we plot A n ,Ai 7 and sinT as functions of Act. As one can see 
there are at least two branches of stable solutions in most cases. In the figures 
the branches are marked with a and b. The third branch that is present in some 
cases is always unstable and we will not be interested in it. The a-branch exists 
for Act < Act^ i.e., in the region where the single- mode fixed-point is unstable, 
as discussed in Sections 4.1 and 4.2. This branch is stable except for a very small 
region near Act = in the case of the I = 1 mode at the first overtone (Fig. 5, left 
panels). We are not sure whether this instability is real or rather an artifact of 
our simplifications. 

The conclusion is that there exists one or two stable stationary solutions for 
all possible values of Act because in the region where a-branch does not exist, 
the single-mode fixed-point is stable. Which one of these solutions, if any, is 



1G 




Act Act 

Fig. 3. Fixed-point solutions for the resonance between the fundamental mode and the nonra- 
dial modes with 1 = 1 (left panels) and 1 = 2 (right panels). The Air-dependence of amplitudes 
(upper panels) and sinT (lower panels) are shown. Solid and dotted lines denote stable and 
unstable solutions, respectively. The range of Act corresponds to A<r max , i.e., to the density of 
the spectra of nonradial modes. The adopted coefficients arc i?=1456, 7 = 0.068, 1 = 0.0811 
for I = 1 and Ft = 1362, 7 = 1.0099, I = 0.0237 for / = 2. 

chosen by the system will be studied by numerical integration of the amplitude 
equations. 

5 Time Dependent Solutions 

In this Section we present time dependent solutions obtained by means of nu- 
merical integration of AEs for the same values of coefficients that were used in 
the stationary solution studies. We focus on the I — 1 cases but the results are 
similar in other cases. 
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Fig. 4. The same as in Fig. 3 but for the nonradial modes with I = 5 (left panels) and I = 6 
(right panels). The adopted coefficients arc _R = 1771, 7 = 4.9738, 7 = 0.0214 for 1 = 5 and 
7? = 2814, 7 = 9.7903, 7 = 0.0265 for 2 = 6. 

5.1 Examples 

Figs. 6 and 7 show the solutions for the resonance between the fundamental 
and I = 1 modes for the selected values of the detuning parameter and different 
initial parameters. In Fig. 6 the detuning parameter of 0.3 is in the range where 
the a-branch exists. The two stationary solutions are AqZh 0.067, 0.013, 
sinr a w 0.142 and A b n 0.059, A\k, 0.022, sinr h « 0.129 (see Fig. 3, left panels). 
Both solutions are stable. In the left panels we show the development of the 
instability beginning with very small pulsation amplitudes. The final (asymp- 
totic) state is an a-branch solution. In the right panels the initial conditions are 
close to a 6-branch solution which is another asymptotic state of the system. 

In Fig. 7 the detuning parameter is 0.65 and, as can be seen from Fig. 3, left 
panels, it is in the range where the a-branch does not exist. Instead we have 
stable single-mode solution Aq = 0.0707, A® = 0. The second stable solution is 
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Fig. 5. The same as in Figs. 3 and 4 but for 
with 1 = 1 (left panels) and 1 = 4 (right panels), 
adopted: R= 1550, 7 = 0.044, 7 = 0.1090 for 1 = 1 



0.76 

0.02 0.04 0.06 0.08 

Act 

the first overtone and the nonradial modes 
The following values of the coefficients were 
and i? = 1748, 7 = 0.1427, 1 = 0.0462 for 1 = 4. 



the 6-branch solution A b x 0.055, A\m 0.026, sinr b = 0.119. The initial condi- 
tions in the left panels are very close to the 6-branch solution and this solution 
is the asymptotic state. In the right panels we start only slightly farther from 
the 6-branch than in the left panels, but the asymptotic state is completely 
different, i.e., it is the single-mode solution. 

As we discussed in Section 4.4 there might be a very narrow range around 
Aer = where a-branch solutions are unstable. We found numerically an am- 
plitude modulated solution with the period of about 60 in dimcnsionless units. 
The range is very narrow and we do not know how robust this finding is. In 
other cases, the asymptotic state is always a stationary solution. Hence we focus 
here only on fixed-point asymptotic solutions. 
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Fig. 6. Time dependent solutions for the same coefficients as used in the left panels of Fig. 3, 
and the detuning parameter 0.3. The initial parameters are Aq = A\ = 0.00001, T = (left 
panels) and A o =0.05, A; = 0.03, T = 0.1 (right panels). 



5.2 Regions of Attraction for the Asymptotic Solutions 

Here we present results of a more systematic survey of the time dependent solu- 
tions. We determine the ranges of the initial conditions (regions of attraction) 
leading to different asymptotic solutions. 

Fig. 8 shows the regions of attraction for the resonance between the funda- 
mental radial mode and the nonradial I — 1 mode. The stationary solutions for 
the adopted coefficients are shown in the left panels of Fig. 3. Fig. 9 shows the 
same solutions in the case of the resonance between the fundamental radial mode 
and the nonradial I = 2 mode. The stationary solutions for the corresponding 
cases are shown in the right panels of Fig. 3. 

We made similar computations for other values of Act and other pairs of 
modes. We found the following properties of the attraction regions. Except for 
the small region near the center of resonance in the case of the first overtone 
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Fig. 7. Time dependent solutions for the same coefficients as used in the left panels of Fig. 3, 
and the detuning parameter 0.65. The initial parameters are Ao = 0.055, A; =0.025, T = 0.12 
(left panels) and A = 0.058, A t = 0.023, T = 0.1 (right panels). 



and 1 = 1 mode, that was discussed in Sections 4.4 and 5.1, all the attractors are 
fixed-point solutions. When o-branch solutions are stable they can be attractors 
simultaneously with a-branch or single-mode stationary solutions. When b- 
branch solutions are unstable then the only asymptotic states are a-branch 
solutions if Act < Act^ or single-mode solutions if Act > Act^. 

From inspection of Figs. 8 and 9 we may conclude that the higher are the 
values of Ai the smaller are the regions of attraction to the o-branch. In other 
words, the farther is the o-branch from the lower (a-branch or single mode) sta- 
tionary solution, the lower is the probability of finding the o-branch pulsation. 
In the three dimensional space of the initial parameters the regions of attraction 
to the o-branch are significantly smaller than those to the a-branch or the single 
mode solution. Let us note also that the star entering the instability strip may 
begin its pulsation only as a single mode pulsator or an a-branch pulsator. Ex- 
cept for the exact resonance, a phase transition to a o-branch pulsation requires 
a finite energy input. Consequently, we will ignore o-branch solutions in the 
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Fig. 8. Regions of attraction for asymptotic solutions for the resonance between the funda- 
mental and the / = 1 modes (see Fig. 3, left panels). Squares denote attraction to the a-branch 
(A<t = 0.3) or to the single-mode fixed-point (Act = 0.6) solution. Dots denote attraction to 
the 6-branch. 
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Fig. 9. Regions of attraction for asymptotic solutions for the resonance between the funda- 
mental and the 1 = 2 modes (see Fig. 3, right panels). Squares and dots denote attraction to 
the a- and fe-branch, respectively. 
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foregoing considerations. 

We see in Figs. 3-5 that if 6-branch solutions are ignored the maximum 
values of the Ai/Ao ratio is at Act = 0. Thus, we may conclude that the strongest 
influence of the resonant nonradial mode excitation on pulsation is at the exact 
resonance. From now on we focus on this special case. 



6 The Role of Nonradial Modes in Constant Am- 
plitude Pulsation 

In this Section we calculate reduction of the radial mode amplitude, Ao/Aq, 
and the relative period change, AP/P, caused by the resonant nonradial mode 
excitation. We also calculate the relative amplitude of the nonradial mode, 
Ai/A , whose presence may perhaps be detected by means of spectroscopy. 



6.1 Analytical Expressions 

As we explained above we assume Aer = 0. Then, if cosT^O, from Eq. (45) we 
immediately get 

(A?) 2 =/(^) 2 , (62) 

or 

f = V7- (63) 

If cosr = then we are not in the a-branch for most of the cases considered. 
The exception is the case in which there is no 6-branch. The example is seen in 
the right panels of Fig. 4. In such a situation the amplitude ratio is determined 
by Eqs.(43) and (44) where we set sinr = l. It happens only in a few cases in 
our large set of resonant mode pairs and we do not consider it in detail. It is 
sufficient to note that the true value of Ai/A is always a bit smaller than \fl 
which we adopt in our calculation. 

In the case of two nonaxisymmetric modes the amplitude of the m = ±l 
modes is given by A± = A\j \/2, which implies 

From Eqs. (43), (44) and Eq. (62) we obtain 

. l+ 7 +2/(g ; -l+ 7 ) 

RSmT s= a l + J(l- 7 ) • (65) 

When sinr s obtained with the above expression is bigger than 1 we should use 
sinr s = 1. This is the case of the absence of the 6-branch discussed above. We 
will not show any numerical results for this case. 
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After inserting solution given by Eq. (65) into Eq. (43) we obtain 
(A s,2 1 1 + 7(1-7) 

(Ao) =a'i+2i(2+my (66) 

which implies 



/ 1 + 7(1-7) ( , 7) 
A° V 1 + 27(2 + 7770' 

For the relative change of the radial mode period caused by nonradial mode 
excitation we get using Eq. (19), 

'^-'^-^Icosr^?) 2 . (68) 



The sign of cosr s is undetermined at Act = 0. In the rest of the a-branch the 
sign of cosT is opposite to the sign of Act. Thus, the frequency shift can be 
either positive or negative with equal probabilities. From Eq. (65) we obtain 



|i?cosr, 



,(1 + 7 + 27(77,-1 +7))' 



\ (i+7(i- 7 )r 



For the dimcnsionless frequency we use 

CT o = — = -5— • (70) 
k Pk 

Finally, with the help of Eqs. (62), (66), (69), (70) we get from Eq. (68) 
IAPJ Pk I 



2ir 1 + 27(1 + 2777, 



x 



— (1 + 7(1 - 7 )) 2 - (1+7 + 27(77, - 1+7)) . (71) 



In a few cases of the lack of the 6-branch the period change is exactly zero. 
Eqs. (67) and (71) allow us to evaluate the relative change of the amplitude 
and period, respectively, for the radial mode caused by a resonant excitation 
of nonradial mode. In the right-hand side we have quantities that may be 
calculated for specified model and mode. The nonradial mode amplitude can 
be calculated with Eq. (63). 



6.2 Application to Evolutionary Models of RR Lyr Stars 

We apply here the formulae of the previous Section to our models of RR Lyr 
stars. Like in Section 4.1 we considered the resonance between the radial and 
the nearest nonradial modes of 1 = 1, 2, 5, 6, for fundamental mode, and of I = 1,4, 
for first overtone. Using the same mode parameters we calculated the nonradial 
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Fig. 10. Finite amplitude consequences of the 1:1 resonance between the radial and the 
selected nonradial modes. Different symbols correspond to different RR Lyr models and 
different modes (see Fig. 1 for more detailed explanation). The top panels show the relative 
amplitudes of the nonradial modes. The middle panels show the ratio of the resonant to the 
nonresonant radial mode amplitudes. The bottom panels show the relative period shift due 
to the resonance. 
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mode amplitudes, Eq. (63), the reduced radial mode amplitudes, Eq. (67) and 
the relative period changes, Eq. (71). Results are given in Fig. 10. The values 
give the maximum effect of the resonance because - we recall - Act = was 
assumed. 

We see in the top panels that the nonradial mode amplitude may be quite 
large. Still detecting the presence of nonradial modes, which would constitute an 
ultimate test for our theory, may not be easy. In the case of low degree nonradial 
modes, the determination of / by means of multicolor photometry is feasible 
(see e.g., Garrido 2000). We should stress, however, that the technique has 
been applied so far only to main sequence pulsators and with only a moderate 
success. For higher degree modes (I > 2) one may contemplate use of line profile 
variations (see e.g., Aerts and Eyer 2000). Again in this case the technique has 
been used only for the main sequence pulsators. 

The amplitude reduction caused by the resonance, which is shown in the 
middle panels, is significant in some cases. Note, in particular, in the right panel 
the relatively large reduction caused by higher degree modes. The implication is 
that the pulsation amplitude must be regarded as a random quantity in certain 
ranges. It is random because the detuning parameter Act which affects the 
amplitude is in fact a random quantity in the [0,ACT max ] range. 

However, in the case of low-degree nonradial modes, the primary source of 
uncertainty of observable amplitude prediction is the orientation of the pulsation 
axis. We will now focus on / = 1 mode and will convert amplitudes shown in 
Fig. 10 to those of light and radial velocity. 

In Appendix B.l we derive formulae for the radial velocity and light curve 
amplitudes. With the help of Eq.(B9) we get the ratio of radial velocity ampli- 
tude with the presence of the nonradial mode, A v , to that without the nonradial 
mode, A v fi. The expression is 

-^ = v/l + e2 + 2e,cos(r/2) (?2) 

where for I = 1 we have 

e v = lMA 1 /A cos& Q . (73) 

Here 9o is the aspect angle, A\/Aq is the calculated amplitude ratio, whose 
maximum values are shown in the top panels of Fig. 10. 
Similarly, we get for the light amplitude ratio 



A 



Mbol 




|cos(?A + r/2) + cos(r/2) 



(74) 



AiWbol,0 

(see Eqs. B12, B19, and B21). For 1 = 1 

e M = l-23Ai/ J 4 cos6 . (75) 

The complex quantity fe 1 ^ is the ratio of the perturbed radiative flux to radial 
displacement at the surface. It is determined from linear nonadiabatic calcula- 
tions. We also denoted 

Po = -i/l + 7cos^ + -^. (76) 
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When we assume the maximum value of the Ai/A ratio to be about 0.3 
the maximum absolute value of e is about 0.4 and 0.37 for the radial velocity 
and the light variations, respectively. The largest effect is for | cos(T/2)| = 1. 
Note that for I = 1 modes we get r s=a or 2ir. Then the rhs in Eq. (72) is 1 ± e v . 
Since e v is a random quantity with values between —0.4 and 0.4, the range of 
the ratio A v /A v ^ ranges from 0.6 to 1.4. 

For the light curve we have to discuss possible values of / and ip- I* 1 our 
linear RR Lyr pulsation models / ranges from about 2 to over 10. The phase 
tp is correlated with / so that it is close to zero for / « 2 and almost tt/2 for 
/ > 10. For the lowest values of these quantities we have po = 2 (see Eq. 76) and 
^MbolMAfboKO = y/l + e 2 M /4 + e M cos(T/2) (see Eq. 74). For |cos(r/2)| = 1 
we obtain Ajtfboi/^Mboi.o = 1 ± €m /2. If / 3> 2, we easily obtain po = 1 and 
^4MboiA4Mboi,o = l± e M, that is almost as wide range as in the case of the radial 
velocity amplitudes. 

At higher Vs, the expected amplitude scatter should be much lower due to 
both the smaller Ai/Aq ratio and the effect of averaging over the stellar disc. 

We conclude that certain scatter in the dependence of the pulsation ampli- 
tudes on stellar parameters should be expected. Let us note that the existence 
of such a dependence underlines empirical methods of absolute magnitude and 
metallicity determination employing light curve data (e.g., Kovacs and Jurcsik 
1996, Jurcsik and Kovacs 1996). 

The period changes shown in the bottom panels of Fig. 10 are indeed rather 
small. Again they should be regarded as random quantities. The expected scat- 
ter is probably inconsequential for the empirical use of period data on RR Lyr 
stars. 

7 Amplitude and Phase Modulation Due to Non- 
radial Mode Excitation 

7.1 Theoretical Prediction 

In Section 4.2 we showed that the stationary solution in the case of the 2:1+1 
resonance is an equally-spaced triplet with the central peak corresponding to 
the radial mode and the two side-peaks corresponding to the nonradial modes. 
We will show here that excitation of the triplet manifests itself, in general, as 
pulsation with amplitude and phase modulation and thus provides a plausible 
explanation of the Blazhko effect. We focus on the I = 1 modes. Modes at 1 = 2 
have significantly lower probability of excitation, as may be seen in Fig. 1, while 
those at / > 2 have significantly lowered amplitudes due to the averaging effect, 
as discussed in Appendix B. We consider first an excitation of the m = ±1 modes 
of the same radial order n. 

In uniformly and slowly rotating stars the frequencies of the m = ±l modes 
are given by 

iu±=uj Q ±cn + Dn 2 , (77) 
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where C (the Ledoux constant) and D depend on the stellar structure. In 
RR Lyr stars contribution to the Ledoux constant arises predominantly in the 
g-mode propagation zone and thus we can use the asymptotic value C = 1/2 for 
1 = 1. Values of D for I = 1 g-modes are very small. As first noted by Buchler et 
al. (1995), the resonant phase-lock effect changes the frequencies in such a way 
that Dil 2 becomes zero. Eq. (77) yields frequencies in the corotating system. In 
the inertial system we have to add — mCl, which accounts for transformation of 
the azimuthal angle. Note that in our notation m = — 1 corresponds to prograde 
modes. Thus, for the nonradial mode frequencies in the inertial system we have 
the following relation 

w±=w TA (78) 

where A = f2/2. Eq. (78) remains valid in the case of depth-dependent rotation, 
but then Q should be understood as the average rotational frequency weighted 
with the Brunt- Vaisala frequency which is strongly peaked in the innermost 
part of the radiative interior. 

In Appendix B.2 we show that when the amplitudes oi 1 = 1, m = ±1 modes 
are equal and small in comparison with that of the radial mode, the radial 
velocity is given by 

V r = A vfl [l + acosAi]cos(o>t-7r/2 + &cosAi) (79) 

(see Eq. B24 and B36), where 

a = e„cos(r/2), (80) 

b = e„sin(r/2). (81) 

The quantity e v is defined similarly as in Eq. (73) but with sin0o instead of 
cos6o, A Vi q is given in Eq. (B6). The maximum value of e v is about 0.4. 
The largest amplitude modulation, factor of 7/3, is obtained for |cos(r/2)| = 1 
and then there is no phase modulation. The opposite case of maximum phase 
modulation and no amplitude modulation corresponds to cos(r/2) = 0. It is of 
interest that for I = 1 we should be close to one of these two extreme cases, since 
we have sinT w 0. 

In a similar manner we obtained (see Appendix B.2 for details) the bolomct- 
ric luminosity variations in the form 

AMboi = AMboi,o [1 + c cos At] cos (tot + po + dcosAt) (82) 

where 

c= — cos(^ + r/2-^o), (83) 
Po 

d=—sm(ij + T/2-^o) (84) 
Po 

(see Eqs. B39, B47, and B48). Quantities p and ip are given by Eqs. (B13) 
and (B14). Explicit expression for p is also given by Eq. (76). For £m we may 
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use rhs of Eq. (75) replacing cos with sin. The lowest values of f,tp are close 
to 2 and 0, respectively. From Eqs. (B13)-(B14) we find po~2, <po~0, which 
yields 

c „ ^cos(r/2), (85) 
d w ^-sin(r/2). (86) 

For the highest values of /(3> 2) we have and we find p ~ 1, ~ 7r /2, 
which yields 

c w e M cos(r/2), (87) 

d w e M sin(r/2). (88) 

This means that we obtain the same conditions for pure amplitude and pure 
phase modulation as in the case of the radial velocity variation. Also the ranges 
of modulation are similar. However, at moderate values of / and ip we may 
have significant both phase and amplitude modulations. 

The 2:1+1 resonance may lead to excitation of certain modes from consec- 
utive multiplets. The two distinct cases involve the following pairs of m values: 
(1,-1), (—1,1). The resonant interaction between a (0,0) pair and the radial 
mode is also possible, but this case is described by much more complicated AEs 
and should be studied separately. Thus we focus here on the m = ±l pairs. 
Then the modulation frequency is 



Aw— Q 



for m's (1,— 1) 



A- 1 AajjtQ for m / s ( _' M) (89) 

where Aw is the frequency difference between centers of multiplets of the radial 
order n and n + 1. In this case we expect unequal amplitudes of the side-peaks 
in consequence of the inequality of mode inertia. In the (—1,1) case the effect 
is systematic. It may be seen in Fig. 2 of DC that lower frequency modes have 
higher inertia typically by a factor of 1.5, implying amplitudes lower by about 
20% (see Eqs. 35, 34). In the (1,-1) case we may have the opposite situation if 
Q>Alj. 

Excitation of nonradial modes provides a natural explanation of basic prop- 
erties of the Blazhko-type modulation. Of course, our simple treatment of non- 
linear effects precludes a detailed comparison with observations. The basic 
characteristics of the modulation are periods and ranges of amplitude variation. 



7.2 Comparison with Observational Data 

We are now in the position to compare quantitative prediction of our model 
with observational data on the Blazhko effect. The numbers to compare are 
modulation periods and ranges of amplitude modulation. 
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The modulation to pulsation period ratio, Pb/Po, ranges from 20 to 1000 
(Kovacs 1993). In our model we have 

ElL = 2 ^L I (90 ) 

where Aw is the frequency difference between centroids of multiplets of radial 
orders n and n+1. We set k = in the case of modes from the same multiplet 
and k — m n in the case of modes from the consecutive multiplets. From Fig. 6 
of DC we find the values of Aw/wo for the 1 = 1 modes in the vicinity of the 
fundamental mode ranging from 4.5 to 7.5 x 10~ 3 . Let us consider first the 
case k 7^0 and f2<C Aw. Then we find 270 < Pb/Pq < 450. Unless our models 
of RR Lyr stars are grossly inadequate, we cannot account for large fraction 
of observed modulation periods. Interpretation of Blazhko stars with short 
modulation periods demands fi>Aw. 

We are not aware of any VJ-otsini determination for any RR Lyr star. Peterson 
et al. (1996) provide only an upper limit of 10 km/s for 8 RRab stars. If we 
assume uniform rotation we may translate this upper limit to the minimum 
value of Pb/Pq- This minimum value is about 30. This is still more than the 
lowest observed modulation period. However, we should remind that f2 does 
not refer to the surface but to the deep interior and it is quite plausible that the 
interior rotates much faster than the surface. We conclude that our prediction 
of modulation periods is not inconsistent with observations. 

A comparison of our prediction regarding ranges of amplitude modulation 
with observational data is much more complicated. Firstly, inadequacy of our 
approximation is certainly far more important in this case; secondly, we have 
much less observational information. The only systematic survey of amplitude 
modulation in RR Lyr stars is that from the MACHO project (Alcock et al. 
2000), but it is restricted to RRc stars in LMC. Let us remind that the Blazhko 
phenomenon is much more common in RRab stars than in RRc. In Alcock's et 
al. (2000) list we find the cases of rather different side peak amplitudes, much 
larger than those implied by our Eq. (64). Also there are cases of side peak 
amplitudes exceeding the central peak amplitude, contrary to our predictions. 
We believe, however, that these cases do not invalidate the very essence of our 
model, but point out to necessity of supplementing it by taking into account 
more interacting modes. Samples of Blazhko RRab stars from the MACHO 
project were analyzed by Kurtz et al. (2000). In a few examples they point out 
different degrees of phase and amplitude modulations. Our model provides a 
natural explanation of the above diversity. 

For a comparison of our simplified theory with observations, the radial ve- 
locity data are more suitable. The reason is that the higher order nonlinear 
effects, which we ignored, are more important in light than in radius varia- 
tions. Unfortunately, the radial velocity data on Blazhko phenomenon are very 
limited. Recently, Chadid et al. (1999) investigated line profile variations in 
RR Lyr itself. They find a triplet associated with the main frequency but it is 
not exactly equally spaced. The spacings are 0.027, 0.020 c/d. We do not know 
whether this departure from equidistancy is significant. The relative side peak 
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amplitudes are about 0.3 which is about the maximum value predicted by our 
model. 

8 Conclusions and Discussion 

We solved simplified amplitude equations describing resonant coupling between 
radial and nonradial modes. The equations cover the case of 1:1 resonance 
between a radial and an axisymmetric (m = 0) nonradial mode as well as the 
case of 2:1+1 resonance between a radial mode and a (to,— m) pair. We showed 
that in the limit of slow rotation and modes of the same multiplet the two 
cases are equivalent. If the modes belong to neighboring multiplets this is not 
true, but the treatment is very similar. The crucial parameter is Aw which 
is a measure of the departure from exact resonance. We showed that if radial 
pulsations are unstable, which requires Aui to be less than certain critical value, 
then there are typically two branches of stationary solutions which merge at 
Alu = 0. The branch corresponding to lower relative amplitude of the nonradial 
mode is almost always stable. The other branch may also be stable in a certain 
range around Alu = 0, but the region of attraction to solutions in the latter 
branch is relatively small. For this reason we use the solutions corresponding 
to Alu — as representing the maximum departure from pure radial pulsation. 

Excitation of an axisymmetric nonradial mode leads to pulsation with con- 
stant amplitude. The presence of a nonradial component affects both period 
and amplitude. The most significant effect is an amplitude change in the case 
of excitation of a dipole mode. The size of the expected amplitude change de- 
pends on Alu and the aspect angle of the pulsation axis. Both must be regarded 
as random quantities. The former, in principle, is determined by stellar param- 
eters, but it is so sensitive that we may specify only its upper limit which is 
(u) nt i — w„+i j ;)/2. We found that for the specified model and pulsation mode 
the uncertainty of the amplitude may reach up to 40%. The consequence is that 
pulsation amplitude and other characteristics of light curve are not uniquely 
determined by mean stellar parameters. This is not a good news for diagnostic 
applications of the light curve data. 

Amplitude modulated pulsation is expected in the case of excitation of non- 
axisymmetric pair of nonradial modes of low degree I. Here again the case of 
I = 1 is most likely and most interesting. In this case we find that the maxi- 
mum range of amplitude modulation is about 40%. The modulation period is 
determined by the characteristic of the deepest layers of the radiative stellar 
interior. Specifically, it is given by 47r/|O rot + k(ou ni i — w„+i,{)|, where fc = 0,±l, 
f2 ro t is the mean rotation rate in those layers, and w n ,! - w n +i,i is the distance 
between the centers of the consecutive multiplets, n at fixed I. Hence, if this is 
the modulation seen as the Blazhko effect then the Blazhko period is a probe of 
deep interior properties of Horizontal Branch stars. The prospect of sounding 
interiors of HB stars is quite interesting. 

The reader must be reminded of simplifications used in this paper. We rely 
on the formalism which assumes that the amplitudes are small. This is not 
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true. Especially our prediction of the light amplitude must be regarded unreli- 
able. Higher order nonlinear effects should be included. Our treatment of mode 
coupling was simplified by adopting a simple scaling of the coefficients and by 
assuming that they are either pure real or pure imaginary. Taking into account 
their complex character may imply an instability of constant amplitude solution 
and imply periodic limit cycle (Moskalik 1986). In this case, the modulation 
period would be determined by the properties of the outermost layers. Validity 
of our formalism is restricted to the cases of 1:1 and 2:1+1 resonances. Fur- 
thermore, couplings which are not covered by our formalism may also play a 
role. These include coupling between the nonradial modes within the multiplet 
as well as the coupling between the radial mode and two axisymmctric nonra- 
dial modes of consecutive multiplcts. We plan to consider these cases in our 
subsequent paper. 

We have realized that there is a considerable paucity of observational data on 
the Blazhko effect. The crucial signature of the type of mode coupling considered 
in this paper is the occurrence of the equidistant triplet in pulsation spectra. 
Such spectra, based on modern photometry of RR Lyr stars, are available only 
for RRc stars and not for RRab stars, for which the Blazhko effect is claimed 
to be much more common. From the point of view of testing our model of 
the Blazhko effect, spectroscopy is more suitable for detecting nonradial modes 
than photometry Further, if such modes are detected we may extract the 
information about mode order and the aspect angle of the pulsation axis. With 
this information we may learn about differential rotation in the interiors of 
Horizontal Branch stars. Spectroscopic studies of Blazhko effect began only 
very recently (Chadid et al. 1999). 
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Appendix 

A Amplitude Equations for Two- and Three- 
Mode Resonance 

In this Appendix wc outline the derivation of amplitude equations for the 1:1 
resonance (the radial and the m = nonradial mode) and the 2:1+1 resonance 
(the radial and the pair ±m of nonradial modes). We briefly sketch how the 
basic form of amplitude equations and the relations between coupling coefficients 
can be obtained. Our main aim is to capture the difference in the coupling 
coefficients for both cases by considering typical third order terms. 

The oscillatory mode is described by the displacement given in Eq. (4). 
The nonlinear equations for the mode amplitudes contain operators acting on 
eigenfunctions and integrated over the volume of a star. For simplicity we will 
be concerned only with the radial component of the displacement. Taking into 
account angular components introduces additional terms in integrals but does 
not change the form of the equations and the relations between coefficients. 
Thus we can write 



tr,j {r,t) = -ctj (t)wj {r,t) + cc, 



(Al) 



where 



w j (r,t) = h j (r)Y l m (0,^ 



iujj t 



(A2) 



Degrees l,m correspond to mode j. Spherical harmonics are expressed by the 
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normalized associated Lcgendrc functions 



Yr(6,<P) = (-ir-}=Pi, m (™s())e im *, (A3) 



and 

i Al+m 

is the associated Legendre function. Normalizations are 

|^, m | 2 df! = l, (A6) 



P,; m (x)di=l. (.47) 

-1 

We will assume the dominant nonlinear terms to be of the third order. These 
terms in equation for the time derivative of the amplitude aj are determined by 
the expression 

T j = Uw j \N\t 1 ,Z 2 ,&), (A8) 
1 i 

where Ij is the mode inertia, N is some operator that is linear with respect 
to each argument and symmetrical with respect to interchanges of arguments 
and 1, 2, 3 correspond to all modes that gives nonvanishing resonant terms. 
The bracket denotes integration over the volume of the star. The operator in 
Eq. (A8) can be presented in the form 

= a 1 a 2 a 3 {h j \N r \h 1 , h 2 ,h 3 ){Y j \N e , !p \Y 1 ,Y 2 ,Y 3 )e i ^ + ^ +Us -^ t . (A9) 

The angular part of the operator may be simplified to the factor J Y*YiY2Y 3 dfl. 

We keep only resonant terms, i.e., terms slowly varying or constant in time. 
This means one of the frequencies cji,W2,W3 has to be taken with minus sign, 
which implies one of the eigenfunctions £1,^2, £3 should be taken complex conju- 
gated. Additionally, integration of the product of spherical harmonics includes 
integration of factor exp[i<^(mi + m 2 + m 3 — m,j)] where one of the numbers 
mi,m2,m 3 is taken with negative sign (complex conjugation of spherical har- 
monic) and for nonvanishing terms the sum of m's in the exponent has to be 
zero. 
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A.l The 1:1 Resonance 

The nonvanishing terms are 

(wo\N\H ,Zo,ti* ) + {w \N\to,$,Zo) + (w \N\%,Z ,t ) = 

= 3(h \N r \h ,h ,h*) [(Y*Y ) 2 dn\a \ 2 a , (AW) 



and similarly 

(w \N\t; ,tit,t;f} + permutations = 

= 6<^o|JV r |/io,^,fer> y|l^| 2 !yi| 2 dn|o,| 2 ao, (All) 

{w Q \N\^,^i,^i} + permutations = 

= 3(ho\N r \h^,hi,hi) / (YjfYfdnc^'-^afa*. (A12) 



All m's are equal to zero. For even / terms involving J F ; 3 do not vanish, but 
we neglect them as was explained in Section 3.1. Symmetrically 

(itfj | iV| + permutations = 

= 3(h l \N r \h l ,h l Xi)J(YiY l *f^l\ai\ 2 a l , (A13) 

(wi\N\£i,£ q ,£q) + permutations = 

= &(h\N r \h u ho,K) J\Yo\ 2 \Yi\ 2 dn\ao\ 2 a u (^14) 

(wi\N\£i ,£ ,£ ) + permutations = 

= 3{h t \N r \ht, h , h ) J Y^Y^dnc-^'-^ala*. (Alb) 

Radial integration may be done only over the envelope. The justification of it is 
that the radial mode cigenfunction is nearly zero in the interior and most terms 
involve it. The term involving only the nonradial mode eigenfunction is the 
saturation term and we are interested only in its real part, as was explained in 
Section 3.1. The real parts of saturation terms are determined by nonadiabacity 
of oscillations which is very weak in the deep interior and it is sufficient to 
restrict integration to nonadiabatic envelope. But in the envelope all radial 
eigenfunctions are nearly the same. It means the brackets on the right-hand 
sides of Eqs. (A10)-(A15) are the same. 
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Let us introduce the quantities 

Joo = 3(h a \N r \h ,ho,h )jY 4 dn^^-(h \N r \h ,h ,h ), (A16) 

Joi = ^ho\N r \h ,h l ,h l )jY 2 Y l 2 dn = ^-(h a \N r \h ,h a ,h a )=J 00 ,(A17) 

Ju = 3(h l \N r \h l ,h l ,h l )jY l 4 dn = 

= ^(ho\N r \h ,ho,h ) J P?(x)dx = 2J oH l , (A18) 

where we used normalization of spherical harmonics, the fact that they are real, 
ho (r) — hi (r) in the envelope, and the easy-to-show property 

1 f 1 

Y 1 Y 2 Y 3 Y 4 dn =—J P 1 P 2 P 3 P 4 dx. 

1 

The quantity Hi is equal J P ; 4 (x)dx. With Eqs. (A16)-(A18) one can write 

-l 

more explicit expressions for the nonlinear terms (A8) 

IoT Q = J o\ao\ 2 a + 2J o\ai\ 2 a + J oa 2 a*c 2[ ^-^ t , (A19) 
IiTi = 2J o\a \ 2 ai + 2J oHi\ai\ 2 ai + J oa 2 atc- 2i ^'-^'> t . (A20) 

With above expressions the amplitude equations are 

^ = Koao + To, (A21) 

da; m . 

= nm+n {A22) 

Relations given by Eq. (13) and the scaling-with-inertia properties (Eqs. 14 and 
15) are obtained in a straightforward manner from Eqs. (A19) and (A20). 



A.2 The 2:1+1 Resonance 

In this case we have nonvanishing resonant terms 



(wo\n\^ + ,c+) 

(wo\N\Zo,Z-,£) 
(w \N\C*,t + ,^) 



3 permutations, 
6 permutations, 
6 permutations, 
6 permutations, 



(w + \N\Z + ,Z + ,Z* + ) 

(w + \N\t + ,£o,$) 
(w + \N\Z+,t-,£) 
(w + \N\£,t ,t ) 



3 permutations, 
6 permutations, 
6 permutations, 
3 permutations, 



37 



and symmetrically for the mode '— '. All properties of the radial eigenfunctions 

are the same as in the previous case. The only m-dependent quantities are 

eigenfrequencies and spherical harmonics. Then quantities given by Eqs. (A16)- 

(A18) are defined analogically, but the product y™y ; ~ m should be used instead 

of Y 2 . Quantity J o is the same as in the m = well as J^l which 

1 

is equal J 00 , too. Finally, J™ = 2J 00 #, m , where HJ n = J P? m {x)dx. Then the 

-l 

nonlinear terms (A8) are 

I T = J ao \a Q \ 2 a + 2J oa (\a + \ 2 + \a^\ 2 )a a + 

+2J oa a-a+a*e^+ +u '-- 2u '^\ (A23) 
I±T± = 2J m \a Q \ 2 a± + 2J 0Q Hl n (\a±\ 2 + 2\a T \ 2 )a ± + 

+J oa a 2 a* T e- l< - UJ + +u >-- 2 ^ t . (A24) 

When we assume equality of nonradial modes inertia, which implies Eq. (24), 
we obtain, similarly as in the previous case, Eqs. (25)-(30). 

B Observed Radial Velocity and Light Curves 

We use here a simple analytical expressions (Dziembowski 1977) which should 
be sufficient for crude estimates presented in this paper. The expression for 
observable radial velocity due to an individual oscillation mode is given by 

(cos 8 ) cos(w/, m i - tt/2 - m$ ) («/ + a H vi) (Bl) 

where l,m are quantum numbers of the mode; u>i ;m is its frequency; i? p h is 
the photospheric stellar radius; A^ m and P^ m are mode amplitude and the 
normalized associated Legendre function, respectively; Go^o determine the 
direction to the observer in the system of coordinates connected with the star; 
ui,vi are certain limb darkening dependent quantities which take into account 
averaging over the stellar disc; and an — G M / (oj 2 R^ h ) . For RR Lyr stars 
this quantity is of the order of 0.1. Similarly we have the expression for the 
bolometric magnitude change 

AMboi = -1.086Ai tm Pi im (<x>sQo)[ficos(ui tm t + ipi-m^o)bi + 

+ cos(w;, TO t-m$ )(26;-Q)] (B2) 

where 6/ , q are again certain limb darkening dependent integrals over the stel- 
lar disc. Quantities fi,ipi are obtained with solving nonadiabatic oscillation 
equations. They allow us to translate the radius changes into the luminosity 
changes. The first term in the square bracket in Eq. (B2) describes emerging flux 
oscillations while the second term describes the radius changes connected with 
oscillatory mode. The values of bi,ci,ui,vi are given by Dziembowski (1977b) for 
Eddington's limb darkening law. For I > 2 they are small due to the averaging 
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effect. The highest values are obviously for 1 = 1 and for radial modes and we 
focus on these modes. Then we have 

P Qi o = y/T/2, Pi,o(cose ) = - N /372cose , 

A,i(c° s ©o) = VW^sinGo. 

The averaging coefficients are u = 0.708, Ui = 0.55, v ~0, «i = 0.45, 60 = 1, 
b\ =0.708, c = 0, and c\ = 1.416. Since aj is small we will neglect terms olrVi. 

B.l The 1:1 Resonance 

Let us consider a superposition of a radial and a m = nonradial resonant 
modes satisfying the double-mode fixed-point solution of AEs. Then their real 
frequencies are exactly equal due to the phase-lock phenomenon. The relative 
phase r is given by Eq. (12) and we see that for phase-locked solution we have 

r = 2(0i iO -0o,o) 
where ifijfi are initial phases in the expression 

9j = -^j-t + (p j fi. 

For individual mode the time dependence is given by cos(u>jt + <pj) = cos(w jt + 
4>j.o), where uij and utj are the linear and nonlinear frequencies, respectively. 
We may choose the initial phase 0o,o to be zero. Then we have <f>jfl = T/2 and 
with the use of Eq. (Bl) for (l,m) = (0,0) and (1,0), we get 

Kad = AoL0 RphUoPo,ocos(uJot-Tr/2) + 

+A l oj l R ph u l P l !0 (cos8 ) cos(uJit -n/2 + T/2) = 
= AowRp h Po,oUo[cos(ut-n/2) + € v cos(u!t-Tr/2 + T/2)], (53) 

where we adopted uj n — uii = ui and introduced 

P (! o(cose ) ui A t 

e v = = — . (±f4) 

P u a A 

Eq. (B3) may be written in the form 

Kad = A t , ! o[cos(wt-7r/2)(l + e t ,cos(r/2))-e t ,sin(r/2)sin(wt-7r/2)] = 

= A Vj0 pcos(ujt — n/2 + Lp) (55) 

where 

A v fi = A u>RphU Pofl (-B6) 
is the radial velocity amplitude for the radial mode, 

pcos^=l + e„cos(r/2) (57) 
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and 

psin<£ = e„sin(r/2). (58) 

From Eqs. (B5), (B7), and (B8) we see that the radial velocity amplitude is 
given by 

A v = A vfi ^fl + + 2e v cos(r/2), (59) 

which is equivalent to Eq. (72). 

With the same choice of phases (fijfi we obtain for the bolometric magnitude 
variations 

AMboi = CA Q P ofi [f cos(ojt + -4;o)bo + cosujt(2bo-c )] + 
+CAiPi,o(cos6 )/i6i cos(ujt + tp 1 +T/2) = 



CA P ,of 



2 

— cosuit + cos(ut + tp) + cm cos(ujt + ip + T/2) 



(BIO) 



where, after Dziembowski (1997b), we assume fi = fo = f and ^pi—tp = tp, for 
our low degree modes. The quantity cm is defined in the same way as e v but with 
6; instead of ui. The radius changes vanish for 1 = 1 modes. The contribution 
to light variation from the radial mode is given by 

AMboi,o = Awboi.o cos(o>i + ipo) (BH) 

where 

A M boi,o = CAoP Qi ofPo, (£12) 



2 

pocosi^o = j+cosip, (B13) 
Posinipo — simp, (£14) 

and 

P 2 o = l + -pr + ycosV-. (£15) 

After adding the contribution from the nonradial mode we obtain 

AMboi = CA P afi f{cosujt[p cos(pQ + e M cos(ip + T/2)}- 
— sinujt [posin<£o + eMsin(^> + r/2)] } = 
= CA P 0fi fpicos(u>t + (pi) = A M boicos(wt + (pi), (B16) 



where 



and 



Picos^i = p cos ipo + e M cos(V' + r/2), (£17) 
pisini^i = p smipo + e M sm(ip + T/2), (BIS) 

A M boi = CA P 0fi fpi. (£19) 
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Eqs. (B13), (B14), (B17), and (B18) yield 

"2 



Pi=Po + e M + 2e M 



cos(V> + r/2)+cos(r/2) 



(B20) 



Using Eqs. (B12), (B19), and (B20) we see that the presence of nonradial mode 
changes the observed amplitude by the factor 




cos(V' + r/2)+cos(r/2) 



(521) 



/ 

This equation is equivalent to Eq. (74). 
B.2 The 2:1+1 Resonance 

We consider here a superposition of a radial mode and a pair of nonradial 
modes. The phase-lock phenomenon in this case causes the three frequencies 
to be equidistant, with the radial mode frequency in the center. The relative 
phase T is given by Eq. (23) and, similarly as in Appendix B.l, for the constant 
amplitude solution, it is given by 

T = 4> m fi + 4>-m,0 - 200,0- 

Again we choose 0o,o = 0. We denote amplitudes and resonant frequencies 
of the radial, the (l,m), and (l,—m) modes as Ao, A_, A + , and ui, u — A, cj + A, 
respectively. Now we have 



rad 



where 



AoLuR ph u Po,o cos(ui - tt/2) + 

+A- (uj - \)Rph,uiPi, m (cos 9 ) cos(wt - At - 7r/2 + (p mfi ) + 
+A + (ui + \)R ph ui Pi, m (cos<d a ) cos(ut + Xt-n/2 + (j>- m ,o) = 
A V fl [cos(ujt — it/2) + e„,_ cos(uit — Xt — it/2 + </>-,o) + 
+e Vi+ cos(o>i + Xt - 7r/2 + 4>+fl)] 

A± u±X ui P (iTO (cos6 ) ui P; !m (cose ) Aj 



(522) 



e v.± = — : ■ — = ~ = : — • (±>Z6) 

Ao lu u Q p 00 u P 00 A 

We made use here of A<w. The amplitude A Vi o is given in Eq. (B6). Keeping 
only linear terms in e Vi ± and making use of certain well-known trigonometric 
relations we get 

V™ d = A vfi [l + a cos(Ai + (3)] cos (ut - tt/2 + 6cos(At + 7)) (524) 

where the modulation coefficients, a, b, and the phase shifts, /3,j, are given by 

acos/3 = e^-Cos^-^ + e^+cos^-^o, (525) 

asin/3 = — e v> _ sin0_ ; o + e„.+ sin0 +! o, (526) 

6sin7 = e v - cos0_, o — £t>,+ cos</> +! o, (527) 

&COS7 = e„ i _sin<^_ ; o + e„ i+ sin0 +i o. (528) 
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For I = 1 and A_ = A + = Ai /\/2 we obtain 

V3 . ui Ai e„, TO 
^± = ^-^00-^=—. (529) 

Inserting it to Eqs. (B25)-(B28) leads to 

acos/3 = ^^-(cos^-^o + cos^+.o) = £ t),mCos(r/2)cos(A/2), (530) 

asin/3 = e„, m cos(r/2) sin(A/2), (531) 
fesin 7 = e„, m sin(r/2)sin(A/2), (532) 
6cos7 = e t ,, m sin(r/2)cos(A/2), (533) 

where T = (f> + ^ + (f>-,o and A = (f> + ^ — Finally we get 

a = e t ,, m cos(r/2), (534) 
b = e„ ;ro sin(r/2), (535) 
= 7 = A/2. (536) 

Since the modulation frequency, A, is much lower than the oscillation frequency, 
u>, the choice of initial modulation phase is unimportant, and we may put A = 
in Eq. (B36) and thus get Eq. (79). 

In a very similar manner we get the following expression for the bolometric 
light variations 



AM{%i = A M boi,o{cosut[cos<p + 



cos At 1,1 i , \ 

H ( cosV'(eM,-Cos0_ i o + eM,+ cos0 +; o)- 



po 



sinV'(eM,-sin^_ i o + eM,+ sin0 +i o)) + 
costp(e M - sin0_, o - £m,+ sin<^ +i0 ) + 
smip(e M ,- cos</>_, - £m,+ cos0 +jO ))] - 



■smojt [sin^o 



cosAi i . , , , n 

H ( sinV>(e M , -COS0_.o + £m.+ cos</> + ) + 

Pq 

+ cosip(e M - sin0_ ;O + £m,+ sin<^ +i0 )) + 

sinAi i . ,, . , . , x 

H — — ( sin^(eM,-Sin0_ iO -eM,+ sin0 +;O )- 

- cosip(e M ,- cos<p-fi - e M ,+ cos^ +i o))] } 

(537) 

where ^-Mboi.o is given by Eq. (B12), po,fo are given by Eqs. (B13), (B14), and 
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For cm, ± Cl we get approximately 

AM b m ol = A Mb oi,o [1 + ccos(Ai + 13')} cos (ut + Lp a + dcos{\t + 7')) (539) 
where 

Po(ccos<A)COS/3' — dsmip cos-f') = 

= e M ,-Cos(V' + <?!)-,o) + eM,+ cos(V' + 0+,o), (540) 

Po(~ccos(po sin/3' + <isin(/?osin7') = 

= eM,-sin(V> + </>-,o)-eM,+ sin(V> + ^+,o), (541) 

Po(csin(/? cos/3' + d cos ipo cos 7') = 

= e M ,-sin(V' + 0-,o) + eM,+ sin(V' + <?!'+,o), (542) 

Po(csini£osin/3' + cicos^osh^') = 

= cm.-Cos^ + 4>- ,o)-eM,+ cos(V' + 0+,o)- (543) 

For A + = A_ =A\I \/2 we have again 

/?'=7' = A/2 (544) 

and then 

ccost^o — rising = — — - cos(V> + T/2), (545) 
Po 

c sin 930 + d cos (p — €M ' m sin(V> + r/2) (546) 

where eM,m = 2em,+ = 2eM,-- The above equations lead to 

c = ^!!i C os(^ + r/2-^o), (547) 

d = ^^sin(^ + r/2-(^o). (548) 
Po 

Again choosing A = in Eq. (B44) we obtain Eq. (B39) to be equivalent to 
Eq. (82) 



